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We discuss the implementation of a multiscale biophysico-chemical model able to cope with the 
main mechanisms underlying cumulative exocytosis in cells. The model is based on a diffusion equa¬ 
tion in the presence of external forces that links calcium signaling and the biochemistry associated 
to the activity of cytoskeletal-based protein motors. This multiscale model offers an excellent quan¬ 
titative spatio-temporal description of the cumulative exocytosis measured by means of fluorescence 
experiments. We also review pre-existing models reported in the literature on calcium waves, protein 
motor activation and dynamics, and intracellular directed transport of vesicles. As an example of 
the proposed model, we analyze the formation of the shield against polyspermy in the early events 
of fertilization in sea urchin eggs. 
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I. INTRODUCTION 

Biological processes frequently involve several mecha¬ 
nisms and sub-processes that link different temporal and 
length scales such as, for example, intracellular molecu¬ 
lar interactions to the scale of cells and beyond to the 
behavior of collectives of cells and even organisms [1]. 
From an applied point of view, biomedical research often 
performs experiments in which these scales are also im¬ 
bricated P]. The quantitative modeling of these systems 
and the processes they perform often requires the exten¬ 
sive use of computational simulations that have to incor¬ 
porate several simulation techniques, each one according 
to the time and space scale associated to the biological 
process in question |3]. 

However, simple theoretical models able to cope with 
these problems are less developed in the literature, al¬ 
though they may be of great importance because they 
can simplify the quantitative description without loss of 
accuracy and, more important, by allowing to give clear 
interpretations of the different mechanisms in terms of 
simple and well understood physicochemical laws. 

In this short review, we aim to describe, link and use 
three theoretical models reported in the literature that 
allow to formulate a theoretical multiscale model able to 
cope, with high precision, with intracellular processes in¬ 
volving exocytosis, that is, with the process by which a 
cell directs the contents of secretory vesicles out of the cell 
membrane and into the extracellular space. We look at 
this particular biological process due to its ubiquity and 
biomedical importance ranging from transcelular trans¬ 
port, insulin, neurotransmitter and enzyme secretions in 
the energetic, nervous and early developmental metabolic 
processes and pathways |1]. In particular, we will study 
the case of the early events of fertilization, that is, the 
formation of the shield against polyspermy in oocytes, 
see Figure |Tj 

Exocytosis is a very broad and complex process that 
involves different sub-processes such as the formation and 
storage of vesicles in different pools inside the cell, the 


calcium signaling process, the active transport, and other 
related processes such as vesicle docking, priming and the 
final fusion of the vesicle in the porosomes of the plas¬ 
matic membraneOlS]. Many of these processes have been 
studied experimentally and, as a consequence, several 
proteins and enzymes have been discovered. However, 
beyond the knowledge of their existence and participa¬ 
tion, there is still a fundamental lack of understanding 
concerning the specific role they play in the secretion 
process. [ 3 . 

Depending on the particular physiological function of 
the exocytosis in each cell type, there are variations in 
the rate of exocytosis, the lag time before it begins, its 
time course, the proportion of vesicles that undergo fu¬ 
sion in response to stimulation, the presence of single 
or multiple granule types, and the nature of the regula¬ 
tion of exocytosis by second messenger pathways [5] . For 
example, in the particular case of fertilization of urchin 
sea eggs, the exocytotic process involves sperm fusion 
with cell membrane, the influx of external calcium and 
the subsequent elevation of calcium concentration due to 
the internal pools. Subsequently, other fundamental pro¬ 
cesses take place like the formation and storage of cortical 
granules, the docking, priming and fusion of the secretory 
vesicles with the plasma membrane, presumably through 
SNARE proteins. Finally, the formation of the vitaline 
envelope and the final endocytosis and recycling of empty 
vesicles are processes that may play an important role in 
regulating the exocytosis [9]. 

The presence of micotubules and actin filaments may 
play opposite roles depending on the cell type. In some 
cases, the excytosis becomes favored by the presence of 
the cytoskeleton, as in the case of serotonin and hyaline 
secretion m, whereas in other cases it may be inhibited 
by the presence of the cytoskeleton that acts as a mechan¬ 
ical barrier nmn]. Presumably, cytoskeleton based se¬ 
cretion depends on the action of protein motors whereas 
the opposite case suggest a diffusion controlled exocyto¬ 
sis m- The distances travelled by the secretory vesicles 
in both cases are very different. In this review, we focus 
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Figure 1: Scheme showing the three main mechanisms associ¬ 
ated to the initial events of oocyte fertilization. Sperm entry 
through the acrosomic reaction liberates the genetic material 
into the oocyte and triggers the local entrance of external 
calcium and the subsequent propagation of a calcium wave 
through the ooplasm. This process activates protein motors, 
that actively transport clusters of cortical vesicles towards the 
plasmatic membrane. Once fused, hyaline is released and ini¬ 
tiates the gradual activation of the vitaline envelope until the 
shield against polyspermy is completed. 

on the first case, that is, in which secretion is mediated 
by cytoskeleton associated protein motors. Therefore, 
for the particular case of the urchin sea fertilization, we 
have considered calcium signaling, vesicle diffusion and 
protein motor activation and its processivity as the rate 
controlling steps in the secretion process. 

An appropriate measure of the cumulative exocytosis 
by release site can be performed with the help of fluores¬ 
cence experiments [MIIllISl- In these experiments, a 
fluorescent dye is dissolved in the culture medium in such 
a way that it is activated when becomes adsorbed in the 
plasma membrane and is illuminated with the appropri¬ 
ate light. During secretion, intracellular vesicles fuse with 
the plasmatic membrane to release its contents. Exposed 
to the culture medium, the membrane of the vesicles ad¬ 
sorbs the dye and increases the local fluorescence, see 
Figure 

In view of this, the theoretical description of cumula¬ 
tive exocytosis is based on an apparently simple model. 
The main question is to know the number of vesicles that 
arrive to a given release site of area A in the membrane 
and fuse with it as a function of time. Therefore, if J(t) 
counts the number of cortical granules per unit area that 
reach the release site at the membrane per unit time, then 
the number n(t) of fused cortical granules as a function 
of time is given by the relation 

n{t) = f 3{dA')dt', (1) 

Jo 



Figure 2: Sample traces of exocytosis in eggs for different 
external calcium concentrations. For each series, the number 
of newly emerged exocytosis disks in each image frame was 
counted. This number plus the numbers from all previous 
frames gave the cumulative number for this time point. Each 
line correspond to different values of external calcium: 0.4, 2 
and lOmM for blue, orange and green points respectively. For 
fitting we use /3 = 1 x w = 223.in equation ([^. 

Other parameters appear in the box inside the text. 

where A is the measure area and d the initial position 
of the vesicle cluster with respect to the plasmatic mem¬ 
brane |15| . The fusion efficiency parameter ^ is the frac¬ 
tion of vesicles that release its content with respect to 
the total number of vesicles arriving at the plasma mem¬ 
brane. Here, for simplicity, we will assume that this 
parameter is the unity. In principle, this simple model 
quantifies the cumulative exocytosis observed in experi¬ 
ments [niin]. 

However, it is clear that the current J (t) must be deter¬ 
mined explicitly in order to make a quantitative progress. 
Determining J(t) is not a trivial task. The mechanism of 
the secretion of the vesicular contents out of the cell is 
a complex process that involves calcium signaling, acti¬ 
vation and action of several cytoskeletal motors through 
distances typically involving few microns. This is an en¬ 
ergy consumption process due mainly to the activity of 
the protein motors and pumps. Therefore, this makes 
very important to consider the detailed description of 
the dynamics and energetics of these entities. 

In fact, in order to obtain 3{t) and perform the pro¬ 
gram outlined, one needs to formulate a multiscale model 
that accounts for the initial events of fertilization. In par¬ 
ticular, we have to reproduce the reaction-diffusion wave 
propagation of internal and external calcium and couple 
it with a directed diffusion model describing the trans¬ 
port of granule clusters containing hyaline to the plasma 
membrane. This second step considers in parallel the 
biochemical kinetics of protein motors, kinesins, myosins 
and dyneins m- As a consequence, our work gives quan¬ 
titative evidence that the long time enzyme secretion of 
cortical vesicles or granules is driven by protein motors. 
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mostly by kinesins. 

The paper is organized as follows. In Section 2 we 
review previous theoretical works on reaction-diffusion 
calcium signaling; protein motor biochemistry, energet¬ 
ics and dynamics; and finally on diffusion in the presence 
of external forces, that allow to propose an unifying view 
of the exocytotic processes from the a physical chem¬ 
istry point of view. In Section 3 we apply the multiscale 
model to describe the slow formation of the shield against 
polyspermy, which occurred when exocyted vesicles lift 
the viteline envelope away. In particular, we join the re¬ 
sults of the pre-existing models on calcium waves, protein 
motors activation and intracellular directed transport of 
vesicles in a cell for a successful quantification of cumula¬ 
tive exocytosis previously reported in experimental data. 
Finally, Section 4 is devoted to present our main conclu¬ 
sions. 


II. REVIEW OF THEORETICAL 
BACKGROUNDS 

As we have already mentioned, the exocytosis of sub¬ 
stances by cells comprise several mechanisms whose in¬ 
terplay is essential. This interplay links several time 
and length scales that should be modeled in a nested 
way. Looking at the events, the entrance of external cal¬ 
cium into the cell through the corresponding ion chan¬ 
nels and its interaction with the intracellular calcium is 
the triggering one of subsequent processes, like activa¬ 
tion of cytoskeletal-based protein motors and their active 
transport of vesicles and granules towards the plasmatic 
membrane along microtubules. Here, we review previous 
theoretical work reported on these three main processes. 

In order to fix ideas, we will discuss the general aspects 
of the three models by considering the phenomenology as¬ 
sociated to the fertilization of sea urchins eggs. This may 
help the reader to create a clearer image of the secretion 
process and the mechanisms involved. 


A. Locally-constrained and directed diffusion 

The mechanisms participating in exocytotic processes 
that we already discussed, calcium wave signaling and 
protein motor action converting chemical energy into me¬ 
chanical one, can be linked by means of a diffusion model 
in the presence of external forces. Because this diffusion 
equation is very general, we will first discuss it in order 
to show where the other mechanisms enter. 

Mass transport inside cells is always affected by ther¬ 
mal agitation that induces a diffusive motion. This dif¬ 
fusion is more or less important depending on the exter¬ 
nal forces acting over the transported bio-particles. In 
addition, the intracellular medium consists of a wide va¬ 
riety of polymers and organelles whose presence causes 


a viscoelastic behavior of the cytosol which, in princi¬ 
ple, should be taken into account El- As a consequence 
of this, when bio-particles undergo different forms of pas¬ 
sive diffusion in the intracellular medium they often show 
anomalous diffusion (subdiffusion) [151 [T^ . 

Subdiffusion has been observed in the passive transport 
of particles in the cellular medium whereas an enhanced 
diffusion has been reported for particles driven by pro¬ 
tein motors |20| . This enhanced diffusion was analyzed 
in terms of a generalized diffusion equation containing 
forces due to the cytoskeleton network and to the pro¬ 
tein motors nnmiisii. However, when protein motors 
carry bio-particles of different sizes, this anomalous dif¬ 
fusion is sometimes weak and normal diffusion can be 
used to describe the motion with an effective diffusion 
coefficient. 

Local confinement and protein motor action are there¬ 
fore the key ingredients to explain vesicle transport in 
cumulative exocytosis [101 1^ . 

Let us consider p(x, t) as the number of transported 
vesicles per unit volume and Dy is the diffusion coefficient 
of vesicles in the cytoplasm. Then, the time evolution of 
p(x, t) is described by the general equation 

|^(x,t)=U„V2p-iv-[pF] (2) 

where F is a vector which represents the total force over 
the transported vesicles and /3 is the friction coefficient of 
the transported vesicles. Here, it has been proposed that 
force F could be a combination of harmonic radial forces 
|101IT51[5T] : 1) a locally constraining force that keeps the 
vesicles fixed to microtubule an 2) a pulling force due to 
the activity of molecular motors. 

Thus, the first force can be modeled as resulting from 
a local harmonic force F = —w^x which is applied simul¬ 
taneously with the force of the motors that move with 
velocity v miniiii. The total force applied on the 
vesicles is an harmonic force displacing in the same di¬ 
rection as the molecular motors with velocity v{t). As a 
consequence of this, it may be modeled by means of the 
relation 


F — TU [x ^niotorij'')] (3) 


where Timotorit) represents the position of the motor as 
a function of time uniiia. Hence, if the protein motor 
velocity is, in general, time dependent, then the force 
takes the form 


F(x,t) = 


(x 


Xo) 



(4) 


where Xo is the initial position of the transported vesi¬ 
cles with respect to the release site. In writing Eq. 
we have supposed that: 1) mobilization of motor does 
not occurred until some critical time to when local Ca^"*' 
has reached some critical value; 2) before its activation. 




4 


kinesin exerts a trapping force —w^(x — xo) over the vesi¬ 
cles, and 3) once activated, the movement and influence 
of the kinesin is only radially outward, being e^. the out¬ 
ward unity vector. Calcium waves and the biochemistry 
of molecular motors, that will be discussed in the follow¬ 
ing subsections, allow to determine an explicit expres¬ 
sion for the average motor velocity v(t) and its activation 
time. In general, we may assume that the process starts 
instantaneously when there exist saturation levels of ATP 
and local Ca^’*' concentration increases over a threshold 
value. This situation can be taken into account through 
a Heaviside function in the form v(t) —>■ 0(t — ti)v{t), 
where ti depends on time and position through Ca^+ 
concentration. 

Inserting Eq. @ in (i) , one can And p using the appro¬ 
priate initial and boundary conditions. The experiments 
suggest that one may use three distinct pools of vesicles 
as initial condition M- one very near membrane, other 
in the center of the cell and one in intermediate position. 
As a boundary condition, one may suppose that vesicles 
reaching the membrane fuse to the plasmatic membrane 
at the same rate at which they arrive. 

From the previous equations, we can know the number 
of cortical granules that reach cell membrane during ex- 
ocytosis. To obtain this quantity, one must calculate the 
current of vesicles J(x, t) implicitly defined by the con¬ 
servation equation pt = V • J. Direct comparison with 
equation Q reveals that the current is given by 

J(x,t) = -D„Vp + /3-VF- (5) 

Since we are interested in cumulative exocytosis, we 
need to sum each vesicle arriving to the membrane as a 
function of time. If J{d,t) counts the density of cortical 
granules reaching the membrane from the initial position 
d of the cluster, then the number n{t) of fused vesicles is 
given by the relation Q. 

As we have mentioned above, we need now to deter¬ 
mine the time at which the motors are activated and the 
time dependence of the molecular motors. In order to 
see how this can be done, in the following subsections we 
will discuss the role and phenomenology of calcium waves 
and signaling, and of the biochemistry and dynamics of 
the molecular motors. 


B. Calcium waves 

Two main models are widely used to explain intracel¬ 
lular oscillations and waves: 1 ) the calcium induced cal¬ 
cium release model (CIRC) and 2) the two-pool model 
(TPM). Both models are based on reaction-diffusion 
equations for the Ca^’*' concentration and one for the 
calcium stored in the internal reservoirs Ca^’*’. The elec¬ 
tion of the model should be based in accordance with cell 
type and the kind of stimulation. By solving these mod¬ 


els with the appropriate boundary condition, one obtains 
the spatio-temporal behavior of the Ca^’*' concentration. 

Inside the cells, after the external stimulation, oscilla¬ 
tions of the calcium concentration produce spatial non¬ 
homogeneities along the cell. These non-homogeneities 
seem to induce a spatio-temporal organization of cell’s 
activity, that is, calcium waves act as a signal starting 
different processes at different times and positions, and 
is a function of the initial stimulation. 

In vertebrates, the Ca^+ is mainly stored in bones, 
from where it may be released to the blood after the ad¬ 
equate hormonal stimulation. Inside many types of cells, 
Ca^“'" it is stored mainly in the endoplasmatic reticulum. 
This allows to use it when necessary. Intracellular Ca^+ 
concentration is lower than the extracellular one, O.lpM 
and IpM, respectively. This implies that the entrance of 
Ca^“'" through the corresponding ion-channels after mem¬ 
brane depolarization is an spontaneous irreversible pro¬ 
cess. This means that no energy is necessary to add in 
order to perform this entrance. After the response of the 
cell to the stimulation is completed, molecular pumps 
regulate again the Ca^’*' concentration inside the cell. 
This one is an energy consuming process which is neces¬ 
sary because Ca^^ may be toxic for high concentrations 
in blood and cells. 

Specifically, the Ca^+ is removed from the cytoplasm 
in two ways: 1 ) pumping it out of the cell and 2 ) pumping 
it inside several organelles like the mitochondria, the en¬ 
doplasmic reticulum (ER) or the sarcoplasmic reticulum 
(SR). Thus, during the initiation and the propagation of 
a calcium wave, the influx of Ca^"*" occurs in two ways: 1 ) 
the entrance of extracellular calcium and 2 ) the release 
of calcium from the internal reservoirs. 

The two internal calcium reservoirs are activated by 
two types of receptors: ryanodine and IP 3 . With the 
attachment of these receptors on the membrane, the ion- 
channels open allowing the release of calcium to the cy¬ 
tosol. Ryanodine is more common in neurons and pitu¬ 
itary cells whereas IP 3 is mainly found in non-muscular 
cells. 

Calcium-Induced Calcium-Release model. The 
calcium-induced calcium release-model (CIRC) was 
originally proposed by Friel in Ref. |23| and it is more 
frequently use in ryanodine dependent processes. 

The CIRC model is based on a balance of Ca^’*' fluxes 
in the cell as it is shown in Figure 2. A single in¬ 
ternal reservoir (ER) exchanges Ca^"*" with the cytosol 
through the fluxes Jl 2 and Jp 2 . The cytosol, in turn ex¬ 
changes calcium with the extracellular environment with 
the fluxes Jli y Jpi- Denoting the calcium Ca^’*' con¬ 
centration in the cytoplasm by c and Cg the calcium con¬ 
centration in the internal reservoir Ca^’*', the model pos¬ 
tulates the following dynamics 
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Figure 3: Flux diagram for the CICR model. In this case, 
a single calcium reservoir is activated due to the increase of 
internal calcium concentration. 


ac 

dt 

I 

II 

+ Jl 2 — Jp2, 

( 6 ) 


dt 

— —Jl 2 + Jp2, 

(7) 

where the fluxes depend on 

the concentrations 

in the 

following way: 

Jli = 

ki{Ce - C), 

Entrance. 


2"^ 

II 

k2C, 

Exclusion. 


Jl2 = 

k3i.Cs - c). 

Release. 

( 8 ) 

Jp 2 = 

k4C, 

Absortion. 



Here, Ce is the external concentration that may be 

assumed as constant due to the large concentration dif¬ 
ference. However, because a linear model do not leads 
to the concentration instabilities responsible for the ap¬ 
pearance of waves and oscillations, it is assumed that 
the reaction rate depends in a non-linear way on the 
cytosolic calcium c in the form 




Ki 


K2C“ 


+ c” 


(9) 


This simple model provides an excellent quantitative de¬ 
scription of experimental cytosolic Ca^'^ oscillations and 
their periods, and predicts the observed flows in each cy¬ 
cle Eg. 

The two-pool model. The two-pool model, originally 
proposed by Goldbeter E3], suggested the existence of 
two calcium reservoirs that are activated in series by dif¬ 
ferent agonists. Specifically, the first reservoir is sensitive 
to IP 3 whereas the second one is sensitive to Ca^+. The 
second pool performs a CICR-exchange process with the 
extracellular medium, see Figure Several peculiarities 
of calcium dynamics for different cellular types can be 
successfully reproduced using this model |24| . 


Outside the cell 


refilling 



Figure 4: Flux diagram for the two-pool model. In this case, 
two different calcium reservoirs sensitive to different agonists 
are activated in series. 


As in the CICR model, a calcium flow balance between 
pools and the exterior is performed. However, in this 
case, the fluxes are given by non-linear relations that 
are determined independently in the experiment. The 
important assumption of the model is that the calcium 
concentration in the first pool (sensitive to IP 3 ) is con¬ 
stant during the whole process. Thus, as in the CICR 
model, only the cytosolic calcium concentration and the 
concentration inside the second pool Cg change with time. 

Hence, assuming also that the IP 3 gives rise to a con¬ 
stant Ca^^ flux r, and that the calcium is pumped to 
the extracellular medium with rate —kc, then the whole 
dynamics is given by 

dc 

— = r - fee-g{c,Cg), (10) 

^=g(c,Cg), (11) 

where g{c,Cs) represents the rate of change of Ca^’*' in 
the second pool that depends on the uptake and release 
flows 

'?(c, Cg) = >fuptake "^release ~ ^/^s, ( 12 ) 

where 

1 - 

'^uptake 

. ^ ( V2CT ] { cP 

release + c™ / V ATf -k c? 

Here, •l^jp^ake which the Ca^+ is pumped 

towards the second pool, that is, by means of an active 
mechanism, and is the rate at which calcium is 

released from this second pool. This is an active feedback 
mechanism for the Ca^+ that is essential for the appear¬ 
ance of calcium oscillations in the cytosol. Finally, kfCs 
is the rate at which the calcium is released to the extra¬ 
cellular medium. The exponents I, m and p and the other 
parameters should be chosen appropriately for each case. 
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Calcium waves and the two-pool model. Frequently, 
Ca^“'" oscillations do not occur in an homogeneous way 
inside the whole cell. On the contrary, these oscillations 
are organized spatially in the form of waves. The velocity 
of these waves is notably similar in many different types 
of cells (5-20/j,to/s). In general, these waves are inde¬ 
pendent of the extracellnlar calcium and they are often 
concentric, plane or even have a spiral shape. The ntil- 
ity of these waves can be attribnted to the intracellular- 
extracellular communication. 

In order to model these waves, it is assumed that the 
cytosolic calcium c diffuses, and therefore it is added a 


diffusion term to the equation (10). Therefore, the re- 


snlting two-pool model for calcium waves takes the form 
dc 

— ^ DcV'^c + r - kc-g{c,Cs), (15) 


^=9(.CCs), (16) 

where Dc is the Ca^’*' diffusion coefficient in the cyto¬ 
plasm |25| . 


C. Biochemistry and dynamics of protein motors 

Kinesins, myosins and dyneis are large families of 
cytoskeleton-based protein motors that participate in the 
exocytosis-endocytosis cycles of cells, among many other 
processes, by transporting organelles and vesicles from 
the inner regions of the cell towards the periphery |26]. 

An important case in which these motors play a key 
role is in the plasmatic membrane resealing of oocytes. 
Here, as we mentioned in the Introduction section, we 
will consider the model case of membrane resealing of 
sea urchin eggs pa HH]. In this case, it has been sng- 
gested that protein motors are involved in the docking 
and delivery of vesicles during exocytosis. Using a con- 
focal microscope, in Ref. m it was observed inhibi¬ 
tion of exocytosis in sea nrchin eggs injected with ki- 
nesin and myosin inhibitors. This observations snpport 
the hypothesis that kinesin and myosin motors mediate 
two sequential transport steps that recruit vesicles to the 
release sites of Ca^^-regulated exocytosis. From this ex¬ 
periment, it was concluded that kinesin leads the slow 
phase of exocytosis delivering vesicles along the micro¬ 
tubules from inner regions of the cell, while myosin leads 
the fast phase, dragging outer vesicles trough filaments 
of F-actin[T3]. It also exists evidence that slow and pro¬ 
longed release of insnlin in /3-cells is also a process medi¬ 
ated by protein motors like kinesins and myosins. Finally, 
it is opportnne to mention that theoretical modeling of¬ 
fered clear snpport to the hypothesis that protein motors 
are responsible for the translocation of vesicles contain¬ 
ing nenrotransmitters in the soma of serotonergic nenrons 

m- 


As proteins, enzymes and molecnlar complexes, the 
processes developed by protein motors nse chemical en¬ 
ergy stored into molecnles such as ATP or GTP, pro¬ 
duced by the mitochondrial system of the cell [T]. This 
is an interesting issne: protein and molecular motors are 
strongly coupled to their environments, from which they 
obtain the "fuel" necessary to perform work during the 
same process they carry ont. That is, molecnlar machines 
have not their inner fuel reservoirs. As a consequence of 
this, their energetics should be described jointly with that 
of the snrrounding medium. 

In vitro experimental stndies determined the depen¬ 
dence of protein motor activity on ATP concentrations 
with well controlled conditions. This approach has ad¬ 
vantages. For instance, it allows for a better analysis of 
the detailed dynamics of the motors irrespective of the 
medinm, but may hide some aspects of its performance 
in in vivo conditions, snch as the transport of proteins, 
RNA, vesicles and even organelles [29] that may be re¬ 
lated, for instance, to several exocytosis-endocytosis pro¬ 
cesses [inill71[Mll3DHS2|- 

Protein motor biochemistry. There are essentially 
two models that represent the motion of protein motors. 
Here, we will focns our attention on the so called hand- 
over-hand model, that was first developed to describe the 
particnlar problem of intracellular transport via kinesin 
motion dne to ATP hydrolysis |33fl36| . 

The biochemical reaction mechanism modeling the ac¬ 
tivity of kinesin motors was originally proposed in Ref. 
(34] and consists on a seqnence of six reactions (see also 
Refs. |371 [55]1. The first reaction describes the forma¬ 
tion of the enzyme-substrate complex MKTa due to the 
capture of an ATP molecule (T) by the microtubule (M) 
kinesin (AT) complex MK with the head a of the kinesin 
attached at position to the microtubule (the head /3 is 
free in a retrograde position): 

MKc + T MKTc,. (17) 

fc_i 

In this scheme, the complex MK plays the role of an 
enzyme that acts over snbstrate T through a catalytic 
reaction |34| . 

The second reaction is controlled by thermal fluctu¬ 
ations that induce conformational changes of the dimer 
constituting the motor stalk. This produces a secondary 
enzyme-snbstrate complex MKT^ that favors a step for¬ 
ward becanse corresponds to an advanced position of the 
/3 head: 


MKTc, ^ MKT'c . (18) 

In general, this step is inflnenced by the load / (cargo 
weight) attached to the kinesin [34]. 

The next reaction describes the attachment of the 
free head domain /3 to the microtnbnle nsing an ADP 
molecnle {D) and forming the complex MKTaDp. This 
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step is considered the slow step of the reaction mechanism 
and therefore the one determining the overall velocity of 
the reaction [33]. It is described by 

MKT'^ % MKr^Dp . (19) 

This reaction occnrs over a high free-energy barrier that 
cannot be snrmounted by thermal flnctnations. After 
this attachment, two reactions occur that prevent the 
detachment of the kinesin from the microtubule 


MKTa^Dp % MKpTa , 

(20) 

MK/sTa MKp{D o P)„ . 

(21) 


The final reaction of the first step corresponds to the 
hydrolysis of ATP at the active site a of kinesin. This 
reaction produces an ADP molecule and an inorganic 
phosphate Pi |3S]. Hydrolysis produces a large amount 
of energy that allows the enzyme to liberate the head a 
and, in this way, surmount the high free-energy barrier 


associated to step (191 


MKp{D o P)„ MKp + D + P, 


( 22 ) 


After this reaction, the cycle is completed and one recov¬ 
ers an initial state MKp one step forward from the initial 
position and with lower free-energy [see Eq. ®1- The 
number of repetitions of this cycle that the molecular 
motor is capable to perform determines its processivity. 

ATP dependence of motor velocities. For our pur¬ 
poses, it is now convenient to consider the dependence 
of protein motors velocity on the ATP concentration in 
the surrounding medium. We first recall that motors 
are activated when cytosolic calcium in its surroundings 
reaches some critical value c^. At this time to, the mo¬ 
tor or the collectivity of motors start their walk along 
the microtubule carrying a cluster of vesicles towards the 
cell membrane. The time of activation depends upon mo¬ 
tor position in the cell with respect to the origin of the 
calcium wave. 

Once the motors have been activated, their average ve¬ 
locity determines the rate at which cortical vesicles, both 
docked to the membrane as far from it, reach an fuse with 
the cell membrane and are exocyted. In Ref. [53] it was 
suggested that this dependence of the velocity on ATP 
follows a Michaelis-Menten (MM) kinetics, and therefore 
the following ATP concentration ([T]) dependence 




'^max [ 7 ^] 

[T] + Km ’ 


(23) 


where Vmax represents the maximum speed of the motors 
and Km is the MM constant of the process. We deter¬ 
mine both constants from experiments reported in Ref. 
|39| for kinesins in sea urchin eggs. In these experiments. 



Figure 5: Kinesin velocity as a function of ATP concen¬ 
tration. Fitting was made using Minimal-Squares with a 
Michaelis-Menten dynamics, eq.(23l. We have found that 
Vmax = 5.25/ims“^ and Km = 14.93^M. 


kinesins adhered to a solid substrate move a filament of 
microtubule for some /rm in presence of several ATP con¬ 
centrations. The velocity of the microtubule is found to 


be [T]-dependent according to Eq. (23l. In the following 
section we will use this dependence to obtain Vmax by 
fitting experimental data with good results. 

ADP inhibition and finite time walks. Several experi¬ 
mental studies showed that ADP may inhibit the proces¬ 
sivity of protein motors, since kinesin, myosin and dyne 
in heads act as multi-substrate enzymes |35l 14011^ . 

More recently, it was probed theoretically in Refs. 
[371155], that taking into account the inhibition by ADP 
in the reaction scheme already discussed, yields a finite 
number of steps of a single walk of a protein motor. In the 
case of kinesins, the average number of processive steps 
predicted was around 60 to 120, in excellent accordance 
with experimental observations. 

Time dependence of protein motors velocity. A con¬ 
sequence of this fact is that the velocity of the protein 
motors is a function of time. This is clear since each step 
uses an ATP molecule that, after hydrolysis produces at 
least one ADP and one Pi molecules. This means: 1) 
The ATP concentration [T] evolves in time following, in 
first approximation, a linear relation 


[T] = -Vmaxt - [r]o, 


(24) 


where [T]o is the initial ATP concentration. 2) Consider¬ 
ing that the autocatalytic production of [ADP] inhibits 
the whole reaction scheme, in particular the formation of 
the MKT enzyme-substrate complex, after many cycles 
of hydrolysis the ADP concentration can be modeled by 
a first order kinetics of the form |37l [35] 

[D]i = [D], + [T]„(l - e-'=-‘). (25) 


After solving the corresponding evolution equations of 
the chemical kinetics, it was shown that the velocity of 







0 


400 



Figure 6: Translational velocity of a kinesin motor as a func¬ 
tion of time given by Eq. ( |26| l with the following values for 
the constants: values: [T]o = ImM, [D]o = 0, = 

900 k^at = 113 = 5,600s“^ and K\j = 

2.24: ^m. Except all values were taken from refer¬ 

ence 1351. Significant values that do not appear explicitly 
in Eq. |26[ ) like Vmax = 11.25/iM s“^, \MK]o = lOOnM and 
K\^ = 28/rM were taken also from reference [35] • Figure 
taken from Ref. |38| . 


a kinesin activated at time to co-n be written in the form 

m 


v{t) 


_ Vmax[T - {t - tp)]_ 

a — -I- [t — (t — to)] ’ 


(26) 


where and a are are parameters that can be measured 
in the experiments and the stopping time is given by 
T = \T]o/vmax- For the cases when ATP concentration 
is low, ADP inhibition controls the dynamics of the ki¬ 
nesin, making its velocity strongly dependent on time, as 
it is shown in Figure]^ However, for usual ATP concen¬ 
trations in living cells, kinesin velocity is practically con¬ 
stant until it stops at a given time that depends on the 
initial ATP concentrations and the catalytic constants. 
In rough terms, the stopping time is around r ^ 70sec, 
a time that can be correlated to that for the duration 
of near pools exocytosis m- Long time exocytosis may 
lasts around few hundreds of seconds. Here, it should 
be stressed that the initial time to is determined by the 
initial position of the molecular motor with respect to 
position at which the calcium wave initiated, and also by 
the calcium wave velocity. 


III. A QUANTITATIVE EXAMPLE: EARLY 
EVENTS OF FERTILIZATION IN SEA URCHIN 
EGGS 

Fertilization in sea urchin eggs has been widely studied 
by biologists of the last century. There is a great amount 
of experimental evidence describing multiple processes 
during activation since the contact of the sperm with the 
jelly: acrosome reaction, calcium wave spreading over 
the cell, exocytosis of cortical granules and formation of 
the fertilization envelope. It is known that a cascade 


of ionic and metabolic changes occur in the sea urchin 
egg after contact with the sperm (i^. The egg response 
to the fertilization stimulus is a rapid depolarization of 
plasma membrane potential establishing the fast block to 
polyspermy at the time that the eggs demonstrate a brief 
global contraction within the first 30 sec after fertilization 

IB]- 

Simultaneously, another signaling process induces the 
production of inositol triphosphate (IP 3 ) and diacyl- 
glicerol. An increase in (IP 3 ) initiates the release of 
Ca^’*' from intracellular stores and produces a calcium 
wave, some seconds after fertilization, that travels across 
the egg form the site of sperm entry to the opposite pole 
|45II46| . This wave of elevated calcium propagates trough 
the entire egg. The maximum of concentration occurs in 
cortical cytoplasm. This wave is preceded by a distinct, 
short transient spike of calcium in the cortex which cor¬ 
relates with the time of sperm-egg binding ED- 

The increase of the intracellular calcium concentration, 
[Ca^“'"]i, stimulates in turn a wave of cortical granule exo¬ 
cytosis, resulting in the elevation of the fertilization enve¬ 
lope and formation of the hyaline layer jUj. It has been 
suggested that this transport of vesicles is carried out by 
protein motors like kinesins and myosins M- Measure¬ 
ments of secreted hyaline, that is, of fused vesicules in 
the plasma membrane has been reported in mi. 

As we have mentioned earlier, kinesin walk is activated 
and regulated by the presence of ATP and intracellular 
calcium HD. In fact, after the local increase of [Ca^+ji, 
kinesins and myosins, powered by the hydrolysis of ATP, 
load the cargo and move along microtubules transport¬ 
ing it from the interior of the cell towards the periphery. 
Using our multiscale model, we show that the calcium 
wave elevates the [Ca^“'']i allowing activation of kinesins 
near the cortex of the cell and mobilizing cortical gran¬ 
ules to the cell membrane where they fuse and spread 
their content to the viteline envelope during the process 
of exocytosis till fertilization envelope is completed. 


A. Results 

Calcium waves. During fertilization, after a period 
of latency which lasts about ~7-40 sec, the eggs of sea 
urchins generate a single ~5-min calcium transient that 
spreads as a wave across the ooplasm HHj. In Lytechinus, 
induced calcium wave is directly preceded by an influx of 
calcium leading to a cortical flash and is followed by a 
few post-fertilization calcium transients m- This wave 
spreads throughout the entire egg with a spherical wave 
front and velocity between 5-10/im/s |49| . 

It is known that external calcium is needed at the on¬ 
set of sea urchin natural fertilization |T|. This calcium is 
available from the sea water in the moment of sperm 
entry during the acrosome reaction EOj. Unfertilized 
sea urchin eggs can release internal calcium stores via 
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Figure 7: Schematic representation of the calcium wave 
spatio-temporal evolution according to the two-pool model. 
Endoplasmatic reticulum is distributed in homogeneous form 
helping to maintain and propagate the calcium wave. Cytoso¬ 
lic calcium diffuses with diffusion coefficient Dc- 


IP 3 or non-IPa -mediated mechanisms [JB]. Following 
insemination, a calcium wave occurs in the presence of 
inhibitors against inositol triphosphate receptors (IP 3 R 
s) or ryanodine receptors (RyR’s) but not when both in¬ 
hibitor types are used, suggesting that the two receptors 
may redundantly participate during fertilization m- 

Therefore, we will assume that the calcium dynamics 
within the fertilized egg can be modeled with the two- 
pool extended model. We will suppose that the endo¬ 
plasmic reticulum (ER) represent the only internal store 
of calcium with the release modulated by both IP 3 and 
Ca^’*'. In addition, we will suppose that the network 
of the ER spreads approximately inside the entire cell 
[5^ . If c and Cs represent concentrations of cytosolic and 
stored respectively, and we assume that only the 

first one is free to diffuse at rate Dc, then the dynamics 
of calcium is dictated by Eqs.(15l and (16). 

We solved Eqs. (15) and (16) using a Finite Element 
Method. As a domain, we used a circle of radius ~50 mi¬ 
crons UH] with a closed border except at the point marked 
by the arrow in Figure |III A| where an initial concentra¬ 
tion slightly higher than the equilibrium concentrations 
simulates the entrance of external calcium due to the 
acrosomic reaction. For the equations (15 16), we fix 
typical values of the chemical parameters and adjust the 
coefficients 0^ and r in such a way that: 1 ) cytosolic cal¬ 
cium spreads as a single wave over the cell, as observed 
in experiments; 2 ) the wave spreads with a circular wave 
front; and 3) the velocity were in accordance with the 
experiments reported in Ref. (51], see also Refs. |44F 
We obtained that the two-pool model reproduces 
remarkably well the spatio-temporal evolution of the cal¬ 
cium wave in the oopplasm. Our results are represented 


Figure 8: The left panel shows the results from our simulation 
of a calcium wave using the Finite Element method. Param¬ 
eters in equations ( |15[ ) and (161 are given by r = 0.0075/iM, 
Vi = 0.45/i M fci = 0.9 mM, Fa = 2.275fj.Ms-\k2 = 2/rM, 
ks = 0.9/iM, k = 0.04s“'^, fc/ = 0.005s“^, Dc = 50pms“^, 
I = 2, m = 2 and p = 4. Cytosolic and stored calcium are set 
up at their equilibrium values c = 0.1666 M and Cr = 2.92M 
in almost all the domain, except by a perturbation in the first 
one of value 1.3/iM in the upper left part of the cell. We have 
used zero flux boundary conditions. Images are taken every 4 
sec. Red and blue correspond to maximum and minimum of 
cytosolic calcium concentration respectively. The right panel 
correspond was taken from Ref. m and shows an experi¬ 
mentally observed calcium wave within a sea urchin egg. 


in the left panel of Figure where the time evolution of 
the calcium wave is compared with the results observed 
in Ref. [ST], right panel. In this figure, at time t = 0, the 
small red semicircle indicates the point at which takes 
place the initial local elevation of calcium concentration, 
that is, where the acrosomic reaction takes places. The 
blue color indicates the low basal calcium concentration 
inside the oocyte before stimulation. Front velocity and 
shape are very similar to the experimental ones. This 
allows to reproduce the formation of the vitaline enve¬ 
lope because the activation of the motors occurs at times 
very similar to the experimental ones. The model does 
not represents well the time scale for the reduction of 
calcium concentration which in the experiments may last 
few minutes. 

Protein motors activation and dynamics. Once hav¬ 
ing the calcium dynamics, then we proceeded to study the 
dynamics of protein motors based on the model discussed 
in Section II.C. One hypothesis of this model is that the 
speed has a Michaelis-Menten dynamics type. In the case 
of kinesin present in the sea urchin egg, this assumption 
can be corroborated through the experiments reported in 
Ref. |39| . In the cited work, a partially purified kinesin 
is adsorbed to a glass coverslip and mixed with micro¬ 
tubules and ATP. The translocating activity of kinesin 
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Figure 9: Possible velocity profiles of molecular motor as a 
function of time . We fix the values of a = 4.8s and /3 = 4.7s 
in equation (26l. The values of kv are 5.6 x 10^s“^ for orange 


a smaller value of this parameter softens engine acceleration 
at the beginning of its walk. The values of Vmax are 5.25, 3.25, 
1.5 and 0.44 micrometers per second from top to bottom. 


to microtubules is observed by contrast microscopy at 
different concentrations of ATP. Experimental data have 
been fitted with Eq. (231 by least squares, corroborating 
that velocity follows a Michaelis-Menten dynamics. 

This result gives the maximum velocity entering Eq. 
(26 \ for the speed of kinesin as a function of time. In 
Figure (|^, we show several possible velocity profiles vary¬ 
ing two different parameters (fc, and Vmax)- While Vmax 
clearly controls the maximum speed reached by a kinesin, 
ky indicates how smoothly the kinesin started his walk. 
In these fittings, we have used that average kinesin walk 
takes about 70 seconds, in accordance with experimental 
results. The molecular motor is rapidly activated after 
calcium increase and starts carrying the transported vesi¬ 
cles to the plasmatic membrane of the oocyte. 


It is convenient to stress that the time at which a mo¬ 
tor starts its walk depends on its distance from the point 
of fertilization and to the cell cortex. Vesicles near the 
entry point of the sperm come quickly to the plasmatic 
membrane. Those initially located far away from the 
stimulation point have a larger latency for the initializa¬ 
tion time ti that depends on the calcium wave velocity 
and shape. The important result represented in Figure 
§ is that, given the ATP ooplasm conditions, the ve¬ 
locity of the protein motors is almost constant within 
the duration of the exocytosis. Note from Figure that 
for large ADP concentrations present, the velocity profile 
as a complicated time dependence, see for instance, the 
dashed-dotted line. 


Directed diffusion and cumulative exocytosis. Finally, 
we quantify the results of cumulative exocytosis reported 
in m by linking the previous theoretical results through 
the diffusion model. In these experiments, the number of 


color 

Dy{y,m^s 

v{fj,ms ^) 


n 

green 

0.6 

1.0 

5 

47 

orange 

0.7 

0.2 

8 

76 

blue 

0.6 

0.05 

8 

27 


Table I: Parameters of the model for the three profiles of Fig¬ 
ure IT] See text for details. 


vesicles that are secreted is counted as a function of time 
over the plasmatic membrane of the egg using confocal 
microscopy. Experiments show that decreasing the ex¬ 
ternal calcium causes a slower rate of exocytosis. These 
results can be understood, at least qualitatively, using 
the diffusion model discussed in Section Hi Al 

In view of the results already discussed for the time de¬ 
pendence of motor velocity, we may simplify the numer¬ 
ical resolution of the diffusion equation by assuming, in 
first approximation, that kinesin moves with a constant 
velocity. Then we solve numerically equation ([^ and 
reproduced the experimental data on cumulative exocy¬ 
tosis by adjusting the parameters (diffusion coefficient, 
friction coefficient, frequency, initial position of the vesi¬ 
cles and motor velocity). As it is shown by the lines in 
Figure |Tj the fits of the cumulative exocytosis with this 
model is excellent, and gives many important information 
on vesicle distribution, motor speeds and the viscoelastic 
properties of the cell that otherwise remain hidden. This 
information is contained in Table I, where the values of 
the diffusion coefficient, the average protein motor veloc¬ 
ities, the initial distance between the clusters of vesicles 
with respect to the plasmatic membrane and the num¬ 
ber of vesicles per cluster are reported. These values 
correlate well with those obtained in other works (cni), 
although in the present case the average velocity is much 
higher, fact possibly related with the biological nature of 
the response. The number of vesicles per cluster is only 
indicative and approximate, since no direct comparison 
with the ultrastructure of the cell allows to compare these 
values. In addition, fluorescence measurements should be 
normalized to a given region and the efficiency of vesicle 
fusion with the membrane was assumed equal to one. 


B. Discussion 

The fits of the calcium wave using the two-pool model 
(with the inclusion of the diffusive term) correctly re¬ 
produce the profile and speed of the calcium wave in its 
initial propagation through oocyte. Notwithstanding, for 
the recovery period in which calcium concentration de¬ 
creases to its basal value we have found (of about 60 
seconds) is well below that found in fluorescence exper¬ 
iments of about 7 minutes |45| . This can be improved 
in our model considering the continuous exchange of cal- 
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cium with the outside and not only at the time of the 
initial application. Furthermore, in our present work, the 
stores of calcium inside the cell are distributed homoge¬ 
neously within the cell, and therefore could no replicate 
the rapid increase of calcium near the cell membrane at 
the moment of fertilization m- We thought this can be 
due to the greater presence of calcium stored in outer 
regions of the egg. However, for qualitative purposes, 
the occurrence of calcium waves founded with the two- 
pool model helps to understand the exchange mechanism 
which produces calcium waves inside the egg. 

We have satisfactory found that engine speed obeys a 
Michaelis-Menten dynamic, corroborating that a kinesin 
can reach a speed of up 5.25fj,ms~^ in in vitro exper¬ 
iments with high concentrations of ATP. However, the 
models presented in this work allow us to think that this 
velocity may be much lower in the fertilized sea urchin 
eggs, where the engine speed when carrying vesicles is 
severely limited not only by the availability of intracellu¬ 
lar ATP [38] , but also by the conditions of friction or drag 
of the intracellular environment where it moves uniini- 
It is therefore not surprising that in our fittings of cumu¬ 
lative exocytosis, the effective velocity we found on our 
adjustments (between 1 and 0.05^ms“^) is substantially 
lower than that of in vitro experiments . For the pur¬ 
poses of this model, we assumed a constant motor speed, 
however, more accurate results may be obtained follow¬ 
ing the time dependence reported in |38| . This may allow 
one to understand that this is only an approximation for 
intracellular processes with a time scale considered com¬ 
parably high to the acceleration and deceleration of ki¬ 
nesin. However, these speed changes can be crucial in 
terms of energy to understand the processivity and stop 
of these motor proteins. 

The diffusion model and drag of the transported vesi¬ 
cles within a cell proposed by Santamari'a-Holek in m 
proves to be extremely successful in quantifying cumu¬ 
lative exocytosis in eggs of sea urchins. However, when 
coupled with the calcium wave patterns, processivity of 
molecular motors, applicability of this model to this bi¬ 
ological system could be expanded to explain and quan¬ 
tify the formation vitaline envelope emerging into the 
egg to prevent polyspermy. This is due to the fact that 
the model perfectly differentiate the start and coupling 
of each process: when the fertilization process starts, it 
emerges a wave of calcium whose concentration is known 
at each position as a function of time |25| : when the cal¬ 
cium concentration reaches locally a critical value at time 
tj, a molecular motor begins its walk dragging outward 
a vesicle |3H]; the rate and the number of vesicles arriv¬ 
ing at each position can be determined by the conditions 
of diffusion and drag of the medium m- Thus, except¬ 
ing the numerical complexity of coupling all the required 
equations, the above model would allow us in the future 
to quantify the spatio-temporal distribution of secreted 
vesicles and understand the formation of the envelope. 


IV. CONCLUSIONS. 

The multi-scale models are very useful when trying to 
describe processes whose causes are described in different 
spatio-temporal scales, fact which occurs in most biolog¬ 
ical processes. In this paper, we have revisited the de¬ 
scription and modeling of three particular physicochemi¬ 
cal mechanisms that are essential in secretion processes: 
a) formation of calcium waves, b) the dynamics of molec¬ 
ular motors based on its biochemistry, and 3) directed 
diffusion of vesicles within a fertilized egg. 

Within this analysis and with the help of previously 
published results in literature, it was found that the ex¬ 
tended two-pool model can be used to understand the 
propagation of calcium waves in the ooplasm of sea urchin 
eggs. Furthermore, the approach we followed in our work 
has the advantage of allowing the integration of specific 
information of other cell types where these waves occur, 
for instance: geometry, exchange flows with the extracel¬ 
lular medium, location of internal pools and intensity of 
initial application. We also reviewed and used a model 
for the processivity of protein motors which is based on 
their biochemistry. This model yields the time dependent 
velocity profiles for the kinesin, protein which plays a key 
role in the translocation of clusters of vesicles in many cell 
types. For this, we reviewed and solved a simple diffusion 
model for vesicle transport that integrates various factors 
affecting the movement of vesicles before fusing with the 
cell membrane, namely diffusion and drag forces due to 
the activity of protein motors. By solving this three com¬ 
ponent model, we have successfully quantified cumulative 
exocytosis previously reported in literature. 

More importantly, in this paper we have described how 
these three schemes can be integrated to quantify with 
accuracy the formation of the viteline envelope. For this, 
in addition to the description of each model, we pro¬ 
posed the relationships between each of them, relation¬ 
ships which are essential to understand this complex pro¬ 
cess. Thus, this work establishes a theoretical foundation 
useful for the understanding and analysis of secretion 
processes, whose difficulty now is merely based on the 
numerical resolution of mathematical equations. 

In future work, we pretend to propose a more detailed 
quantitative model which will consider aspects related 
to the kinetics of docking and priming of the vesicles in 
the porosomes EllSI, as well as the detailed mechanical 
description of swelling of the vesicles and the expulsion 
of its content. 
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